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Abstract 



Zero-point fluctuations in quantum fields give rise to observable forces be- 
tween material bodies, the so-called Casimir forces. In these lectures I present 
the theory of the Casimir effect, primarily formulated in terms of Green's func- 
tions. There is an intimate relation between the Casimir effect and van der 
Waals forces. Applications to conductors and dielectric bodies of various shapes 
will be given for the cases of scalar, electromagnetic, and fermionic fields. The 
dimensional dependence of the effect will be described. Finally, we ask the ques- 
tion: Is there a connection between the Casimir effect and the phenomenon of 
sonoluminescence? 

1 Introduction 

The Casimir effect is a manifestation of macroscopic quantum field theory. Its origin 
is intimately tied to the van der Waals forces between neutral molecules. In 1873 van 
der Waals [|IJ introduced these weak intermolecular forces to explain derivations from 
the ideal gas laws, but their physical basis did not begin to emerge until 1930 when 
London [T2J showed how the leading behavior of this force could be understood from 
quantum mechanics. That is, the Hamiltonian describing the interaction between two 
dipoles separated by a distance R is 

di • d2 — 3dx • n d 2 ■ n 
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where n is the relative direction of the two dipoles. Because (di) = (d2) = 0, this 
Hamiltonian has vanishing expectation value, (H int ) = 0; however, in second order of 
perturbation theory this operator appears squared, so there is an effective interaction 
potential between the dipoles, 

V eS ~ -L (1.2) 

This is a short distance electrostatic effect. The next step in the saga occurred in 1948 
when Casimir and Polder H introduced retardation, resulting in a potential which 
for large distances falls off faster, 

Kff ~ ^. (1.3) 

This result can be understood by dimensional considerations. For sufficiently weak 
fields, the polarization of a molecule should be proportional to the electric field acting 
on the molecule, 

d = «E, (1.4) 

where a is the polarizability of the molecule. Each atom produces a dipole field, so the 
two atoms polarize each other. In terms of a unit of length L, we have dimensionally 

[d] = eL, [E] = eL' 2 , so [a] = L 3 . (1.5) 

Thus we have, at zero temperature, 

«i«2 he 

V eS ~ - sr -. (1.6) 

[It may be useful to remember the conversion factor he ~ 2 x 1CT 5 eV cm.] The high 
temperature limit is classical, 

V eS ~^kT, T^oo. (1.7) 

Here we have thought of the interactions from an action-at-a-distance point of 
view. From this viewpoint, the Casimir effect is simply a macroscopic manifestation 
of the sum of a multitude of molecular van der Waals interactions. However, shortly 
after the original Casimir and Polder paper, Bohr suggested to Casimir that zero-point 
energy was responsible for the intermolecular force, and indeed that is an equivalent 
viewpoint [|J . One can shift the emphasis from action at a distance between molecules 
to local action of fields. 

The connection between the sum of the zero-point energies of the modes and the 
vacuum expectation value of the field energy may be easily given. Let us regulate the 
former with an oscillating exponential: 

\ J2 hu a e-^ = | r ^ e -^ J2 2 \ ■ , r -> 0, (1.8) 
2 ' 2 J_ 00 2m jL ^ / uj l a — uj l — le 
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where a labels the modes, and the contour of integration in the second form may be 
closed in the lower half plane. For simplicity of notation let us suppose we are dealing 
with massless scalar modes, for which the eigenf unctions and eigenvalues satisfy 

-VVo = "fya- (1.9) 

Because these are presumed normalized, we may write the second form in Eq. (|1.8| ) 

as 



u; 2 — ze 



(rfx)9 o 9' o (0(x)0(x'))L„_, (1.10) 



where the Green's function G(x, £; x', i') satisfies 

-V 2 + ^ G(x, f ; x', = *(x - x')5(t - 0, (1.H) 

and is related to the vacuum expectation value of the time-ordered product of fields 
according to 

G(x,i;x',0 = i(T0(x,t)0(x',t')>. (1-12) 
For a massless scalar field, the canonical energy-momentum tensor is 

= <9><9> - ^<9><9 A 0. (1.13) 

The second term involving the Lagrangian in Eq. (|1.13|) may be easily shown not 
to contribute when integrated over all space, by virtue of the equation of motion, 
—d 2 (f) = outside the sources, so we have the result, identifying the zero-point energy 
with the vacuum expectation value of the field energy, 

^X>.= /(dx)(T 00 (x)>. (1.14) 

a 

In the vacuum this is divergent and meaningless. What is observable is the change 
in the zero-point energy when matter is introduced. In this way we can calculate the 
Casimir forces. 

Alternatively, we can calculate the stress on the material bodies. Figure [1] shows 
the original geometry considered by Casimir, where he calculated the quantum fluc- 
tuation force between parallel, perfectly conducting plates H. The force per unit 
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Figure 1: Geometry of parallel conducting plates. 

area / on one of the plates is given in terms of the normal-normal component of the 
stress tensor, 



/ 



dx dy(T zz ), 



1.15) 



where the integral is over the surface area of the plate in question. For electromagnetic 
fields, the relevant stress tensor component is 



T, 



-(Hi- Hi + El- El). 



We impose classical boundary conditions on the surfaces, 



H, 



0, E 



0. 



;i.i6) 



;i.i7) 



and the calculation of the vacuum expectation value of the field components reduces 
to finding the classical TE and TM Green's functions. In general, one further has to 
subtract off the stress that the formalism would give if the plates were not present, 
the so-called volume stress, and then the result of a simple calculation, which is given 
below, is 



f = J dxdy [T zz - T zz (vo\)] 



rr 



240(1* 



he. 



;i.i8) 



an attractive force. 

The dependence on a is, of course, completely determined by dimensional consid- 
erations. Numerically, the result is quite small, 



./ 



-8.11MeVfma" 



-1.30 x 10" 27 Nm'a 



;i.l9) 



and will be overwhelmed by electrostatic repulsion between the plates if each plate 
has an excess electron density n greater than 1/a 2 , from which it is clear that the 
experiment must be performed at the /im level. Nevertheless, it has been directly 
measured to an accuracy of several percent || [7], [|, |, |K| [TT], [12], [13], [L4|, [15|] . (The cited 
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measurements include insulators as well as conducting surfaces; the corresponding 
theory will be given in Section 3.) Until recently, the most convincing experimental 
evidence comes from the study of thin helium films jnj|; there the corresponding 
Lifshitz theory fT/] , 18] has been confirmed over nearly 5 orders of magnitude in 
the van der Waals potential (nearly two orders of magnitude in distance). Quite 
recently, the Casimir effect between conductors has been confirmed to the 5% level 
by Lamoreaux |19|, and to 1% by Mohideen and Roy pO . 



2 Casimir Force Between Parallel Plates 
2.1 Dimensional Regular izat ion 

We begin by presenting a simple, "modern," derivation of the Casimir effect in its 
original context, the electromagnetic force between parallel, uncharged, perfectly con- 
ducting plates. No attempt at rigor will be given, for the same formula? will be derived 
by a consistent Green's function technique in the following subsection. Nevertheless, 
the procedure illustrated here correctly produces the finite, observable force starting 
from a divergent formal expression, without any explicit subtractions, and is therefore 
of great utility in practice. 

For simplicity we consider a massless scalar field <fi confined between two parallel 
plates separated by a distance a. (See Fig. [1].) Assume the field satisfies Dirichlet 
boundary conditions on the plates, that is 

0(0) = <P(a) = 0. (2.1) 

The Casimir force between the plates results from the zero-point energy per unit 
transverse area 



2 ^ J (2^)2 V 



where we have set U — c — 1, and introduced normal modes labeled by the positive 
integer n and the transverse momentum k. 

To evaluate Eq. ( j2.2| ) we employ dimensional regularization. That is, we let the 
transverse dimension be d, which we will subsequently treat as a continuous, complex 
variable. It is also convenient to employ the Schwinger proper-time representation 
for the square root: 

where we have used the Euler representation for the gamma function. We next carry 
out the Gaussian integration over k: 



1 1 ^ f°° dt 



u 



4V5F(47r)"/2^y t 
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Finally, we again use the Euler representation, and carry out the sum over n by use 
of the definition of the Riemann zeta function: 

..-''(f'rf-^d-i-D. (2.5) 

When d is an odd integer, this expression is indeterminate, but we can use the reflec- 
tion property 

r (|) a^- z/2 = r (^) C(i - ^ ( ^ 1)/2 (2.6) 



to rewrite fl2T50 as 



"=-^7^^ r ( 1 + C(2 + <() - P7) 

We note that analytic continuation in d is involved here: Eq. ( |2.4j ) is only valid if 
Re d < — 1 and the subsequent definition of the zeta function is only valid if Re d < —2. 
In the physical applications, d is a positive integer. 

We evaluate this general result (|2.7|) at d — 2. This gives for the energy per unit 
area in the transverse direction 

7T 2 1 

^"1440^ ^ 

where we have recalled that £(4) = 7r 4 /90. The force per unit area between the plates 
is obtained by taking the negative derivative of u with respect to a: 

! - = -l u = "55?- < 2 ' 9) 

The above result (|2.9| ) represents the Casimir force due to a scalar field. It is 
tempting (and, in this case, is correct) to suppose that to obtain the force due to 
electromagnetic field fluctuations between parallel conducting plates, we simply mul- 
tiply by a factor of 2 to account for the two polarization states of the photon. Doing 
so reproduces the classic result of Casimir 0: 

TT 2 1 

f em = (2.10) 

J 240 a 4 V ; 

A rigorous derivation of this result will be given in Sec. |2.5| . 
2.2 Scalar Green's Function 



We now rederive the result of Sec. |2J] by a physical and rigorous Green's function 
approach. The equation of motion of a massless scalar field <j) produced by a source 
K is 

-d 2 (f) = K, (2.11) 
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from which we deduce the equation satisfied by the corresponding Green's function 

-d 2 G(x,x) = 8{x-x'). (2.12) 

For the geometry shown in Fig. [I], we introduce a reduced Green's function g(z,z') 
according to the Fourier transformation 



G(x, x) 



d k ik-(x-x') f -iu(t-f) 

(2ir) d J 2tt 



(2.13) 



where we have suppressed the dependence of g on k and u>, and have allowed z on the 
right hand side to represent the coordinate perpendicular to the plates. The reduced 
Green's function satisfies 



dz 2 



A 2 ) g(z, z) = 6(z-z'), 



(2.14) 



where A 2 = u 2 — k 2 . Equation ( |2.14j ) is to be solved subject to the boundary conditions 
(|1|), or 

g(0 : z')=g(a,z') = 0. (2.15) 
We solve (|2.14j) by the standard discontinuity method. The form of the solution 



is 



AsmXz, < z < z' < a, 

B sin \(z — a), a > z > z' > 0, 



(2.16) 



which makes use of the boundary condition (|2.15|) on the plates. According to 
Eq. ( |2.14| ), g is continuous at z = z', but its derivative has a discontinuity: 



A sin \z' — B sin X(z' — a) = 0, 
AX cos Xz' — BX cos X(z' — a) = 1. 

The solution to this system of equations is 

1 sin X{z' — a) 



A 



B = - 



X sin Xa 
1 sin Xz' 



X sin Xa 

which implies that the reduced Green's function is 

, „ 1 



g\z,z 



X sin Xa 

where z > (z<) is the greater (lesser) of z and z' 



sin Az< sin X(z > — a) 



(2.17a) 
(2.17b) 



(2.18a) 
(2.18b) 

(2.19) 
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From knowledge of the Green's function we can calculate the force on the bounding 
surfaces from the energy-momentum or stress tensor. For a scalar field, the stress 
tensorQ is given by Eq. ( |1.13| ). What we require is the vacuum expectation value of 
Tju, which can be obtained from the Green's function according to Eq. ( 1.12 ), or 



{(j){x)(j){x')) = -G(x,x') 



(2.20) 



a time-ordered product being understood in the vacuum expectation value. By virtue 
of the boundary condition (|2.1| ) we compute the normal-normal component of the 
Fourier transform of the stress tensor on the boundaries to be 



(T zz ) = ^d z d z ,g{z, z') 



2— >z'=0,a 



-A cot Aa. 
2 



(2.21) 



We now must integrate on the transverse momentum and the frequency to get the 
force per unit area. The latter integral is best done by performing a complex frequency 
rotation, 



to 



Thus, the force per unit area is given by 



ip. 



f 



d d k 



(2vr) d J 2tt 



pcoth pa. 



(2.22) 



^2.23) 



This integral does not exist. 

What we do now is regard the right boundary at z = a, for example, to be a 
perfect conductor of infinitesimal thickness, and consider the flux of momentum to 
the right of that surface. To do this, we find the Green's function which vanishes at 
z — a, and has outgoing boundary conditions as z 
like that which led to Eq. (|2.19| ) yields for z, z' > a 



oo, ~ e thz . A calculation just 



g{z,z') = - [sm\( Z< - a)e iX{z >- a) ] . 



(2.24) 



(This obviously has the correct discontinuity — ■§^g\ z z ,~^ ) = 1 and satisfies the appro- 
priate boundary conditions.) The corresponding normal-normal component of the 
stress tensor at z = a is 



(T zz ) 



—d z d z ,g(z,z') 



z=z'=a 



A 
2' 



(2.25) 



1 The ambiguity in defining the stress tensor has no effect. We can add to T^ v an arbitrary 
multiple of {d^d v — g^ v d 2 )(j) 2 21, ^2|. But the zz component of this tensor on the surface vanishes 
by virtue of Eq. (2.1). 



8 



So, from the discontinuity in T zz , that is, the difference between fl2.21| ) and (|2.25|) , we 
find the force per unit area on the conducting surface: 



/ = -; 



d d k f dc 



(27T) 



2ty 



p(coth pa — 1). 



(2.26) 



We evaluate this integral using polar coordinates: 



/ 



A 



d+l 



(27T 



, d+l 



p d dp 



p 



g2pa I ' 



(2.27) 



Here A n is the surface area of a unit sphere in n dimensions, which is most easily 
found by integrating the multiple Gaussian integral 



d n xe~ x2 = ir n/2 



in polar coordinates. The result is 



A, 



2TX n ' 2 



n r(n/2)' 

When we substitute this into (|2.27|) , employ the integral 



dy 



ey - 1 



r(s)((s) 



and use the identity 



T(2z) = (2rr)- 1/2 2 2z - 1/2 T (z) T + 



we find for the force per unit transverse area 



f = -(d + l)2- d - 2 ir- d/2 - 



_ 1 T(l + l)((d + 2) 



,d+2 



Evidently, Eq. ( |2.32|) is the negative derivative of the Casimir energy 
to the separation between the plates: 



du 

da ' 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.32) 
with respect 

(2.33) 



this result has now been obtained by a completely well-defined approach. The force 
per unit area, Eq. (p.32[) , is plotted in Fig. ^ where a — > 2a and d = D — 1. This 
result was first derived by Ambj0rn and Wolfram . 
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-10.0 



-5.0 -3.0 -1.0 1.0 3.0 5.0 

D 

Figure 2: A plot of the Casimir force per unit area / in Eq. ( |2.32| ) for — 5 < D < 5 
for the case of a slab geometry (two parallel plates) . Here D — d + 1 . 



We can also derive the same result by computing the energy from the energy- 
momentum tensor.0 The relevant component is^ 

Too = l(do<pdo<P + d 1 <pd 1 <p + «9 2 0«9 2 + <9 3 0<9 3 0), (2.34) 
so when the vacuum expectation value is taken, we find from Eq. ( |2.19| ) 

(Too) = ; r~ [(^ 2 + k 2 ) sin Xz sin X(z — a) 

2i\ sin Xa 

+ X 2 cos Xz cosX(z — a)] 
'uj 2 cos Xa — k 2 cos X(2z — a)]. (2.35) 



2iX sin Xa 



We now must integrate this over z to find the energy per area between the plates. 
Integration of the second term in Eq. ( [2.35 ) gives a constant, independent of a, which 
will not contribute to the force. The first term gives 



ijj 2 a 



dx (Too) — ~~2 7 \ co ^ (2.36) 



2 Again, the ambiguity in the stress tensor is without effect, because the extra term here is V 2 ^ 2 , 
which upon integration over space becomes a vanishing surface integral. 

3 As noted after Eq. ( 1.13 ), we would get the same integrated energy if we dropped the second, 
Lagrangian, term in Too there, that is, used Too = do4>da(f>. 
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As above, we now integrate over u and k, after we perform the complex frequency 
rotation. We obtain 



u— — / 7—7-7 / — — cothpa. (2.37) 
2 J {2Ti) d J 2rr p H K ' 

If we introduce polar coordinates so that Q = pcosd, we see that this differs from 
Eq. ( p.23|) by the factor of a(cos 2 #). Here 

, „ 4 C cos 2 9 sin**" 1 9d9 1 

cos 2 9) = J ° r . d _ lQ ., = j—, 2.38 
J sin 9d9 d + 1 

which uses the integral 

[\m d - l 9d9= [\l- x 2 Y d - 2 )/ 2 dx = 2 d - 11 -^-. (2.39) 
Jo J-i r(a) 



Thus, we again recover Eq. 

For the sake of completeness, we note that it is also possible to use the eigenfunc- 
tion expansion for the reduced Green's function. That expansion is 



2 sm(nnz/a) sm(nrrz' /a 

n=l 

When we insert this into the stress tensor we encounter 



. Z v—v sm{nirz/a) simnirz /a) .„ , . 

a ^—f n z 7T z /a z — A z 



„„ , 2 ^ n 2 ir 2 /a 2 , n „ 1N 

n=l 

We subtract and add A 2 to the numerator of this divergent sum, and omit the diver- 
gent part, which is simply a constant in p. Such terms correspond to 5 functions in 
space and time (contact terms), and should be omitted, since we are considering the 
limit as the space-time points coincide. We evaluate the resulting finite sum by use 
of the following expression for the cotangent: 

12x^1 , 
COt 7T.T = — + — > — -. 2.42 

k=l 

So in place of Eq. Q2.21 ) we obtain 

(T zz ) = ^ A ( cotAa "^) > ( 2 - 43 ) 
which agrees with Eq. (|2.21|) apart from an additional contact term. 
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2.3 Massive Scalar 



It is easy to modify the discussion of Sec. |27| to include a mass \i for the scalar field 



The reduced Green's function now satisfies the equation 

Q2 



\ 2 + ^ 2 )g(z,z>) = 6(z-z f ), (2.44) 



dz 2 

instead of Eq. ( |2.14j ), so the reduced Green's function between the plates is 

g(z,z>) = _l^., < sin^ > -a) 
k sin Ka 

where 

k 2 = X 2 — /x 2 . (2.46) 
The calculation proceeds just as in Sec. |2.2| , and we find, in place of Eq. ( [2. 27j ) 



f = --^±-r p « dp yf£±? ■ (2.47) 



When we substitute the value of A^+i given by Eq. ( [2.29|) , and introduce a dimen- 
sionless integration variable, we find for the force per unit area 

1 f°° 2 (x 2 -4/i 2 a 2 )( rf - 1 )/ 2 

J 2 2 («)#i)/2r(f)a^t e'-l ' [ ] 

For d = 2 this function is plotted in Fig. The massive scalar was first treated 
by Hays in two dimensions ||24j| , and again done in general number of dimensions by 



Ambj0rn and Wolfram |23 



2.4 Finite Temperature 

We next turn to a consideration of the Casimir effect at nonzero temperature. In 
this case, fluctuations arise not only from quantum mechanics but from thermal 
effects. In fact, as we will shortly see, the high-temperature limit is a purely classical 
phenomenon. 

Formally, we can easily obtain the expression for the Casimir force between parallel 
plates at nonzero temperature. In Eq. (|2.26|) we replace the imaginary frequency by 

C - t = (2.49) 

where (3 = 1/kT, T being the absolute temperature. Correspondingly, the frequency 
integral is replaced by a sum on the integer n: 



2tt (3 ^ 
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2.0 3.0 
Ha 



Figure 3: Scalar Casimir force between parallel plates as a function of mass for d = 2. 



Thus, Eq. Q2.26Q is replaced by 



d d k 



2(3 J (2tt) 



E 



2p r , 



n=— oo 



(2.51) 



where p n = \Jk 2 + (2irn/f3) 2 . 

We first consider the high-temperature limit. When T — >• oo {(3 —>■ 0), apart 
from exponentially small corrections, the contribution comes from the n = term in 
the sum in Eq. (|2.51|) . That integral is easily worked out in polar coordinates using 
Eqs. ( ggg ) and ( PD| ). The result is 



./ 



T^oo 



—kT- 



d 



i ■[*±±)c(d + i) 



(2^a) d + l ~ V 2 

In particular, for two and three dimensions, d — 1 and d = 2, respectively, 



d=l: / 
rf = 2: / 



-kT 
-kT 



71 



24a 2 
C(3) 
8vra 3 



(2.52) 

(2.53a) 
(2.53b) 



This high-temperature limit should be classical. Indeed, we can derive this same 
result from the classical limit of statistical mechanics. The Helmholtz free energy for 
massless bosons is 



-kTlnZ, \nZ 



-ppi 



(2.54) 
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from which the pressure on the plates can be obtained by differentiation: 

dF 



P = ~ 



dV 



^2.55) 



We make the momentum-space sum explicit for our d + 1 spatial geometry: 



F = kTV 



J^y^l E ln(l-exp(-^ 2 + n 2 7r 2 /a 2 ))- (2-56) 



Now, for high temperature, (3 —* 0, we expand the exponential, and keep the first 
order term in j3. We can write the result as 



F = VkT±-± 



d d k 



2a ds J (27r 



t \r ( 



n 2 7r 2 



(2.57) 



s=0 



where we have used the identity 



ds 



(2.58) 



s=0 



This trick allows us to proceed as in Sec. |2.1| . After the k integration is done, the s 
derivative acts only on l/r(— s): 



d 



ds r(— s 

so we easily find the result from Eq. (ET 



(2.59) 



s=0 



(2.60) 



The pressure, the force per unit area on the plates, is obtained by applying the 
following differential operator to the free energy: 



_d_ 

"dV 



]_d_ 



(2.61) 



where V = Aa, A being the <i-dimensional area of the plates. The result of this 
operation coincides with Eq. (|2.52|) . 

The low-temperature limit (T — >• or (3 — >• 00) is more complicated because f T is 
not analytic at T = 0. The most convenient way to proceed is to resum the series in 
Eq. (|2.5ip by means of the Poisson sum formula, which says that if c(a) is the Fourier 
transform of b(x), 



c(«) = 7^ / b(^)e 



dx, 



(2.62) 
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then the following identity holds: 

OO CO 

b(n) = 2tt c(2vm). (2.63) 



Here, we take 



6 M = /(Skfe' ' 2 ' 64 ' 

Introducing polar coordinates for k, changing from k to the dimensionless integration 
variable z = 2ap x , and interchanging the order of x and z integration, we find for the 
Fourier transform 

c («) = 2W^i dxcosax\z -{—) j ■ (2-65) 

The x integral in Eq. ( |2.65| ) is easily expressed in terms of a Bessel function: It is 
^ fl du cos ( ^u) (1 - u 2 f~V 



Ana ./n V Ana 

(d-3)/2 



a (l-d)/2 z (d-l)/2 



ina 



We thus encounter the z integral 



;°° dz 2 ( d + 3 )/ 2 

J ( s ) = / eZ _ 1 WK ( 2 - 67 ) 



where s = a/3/4.7ra. 

The zero-temperature limit comes entirely from a = 0: 





So we require the small-x behavior of the Bessel function 



f=° = -£c(0). (2.68) 



/ Z \ (<i-l)/2 1 



whence using Eq. (|2.30|) 

,27 v^F 2 



\ (d-i)/2 9 d+1 d + 1 /r/ 
1(a) (? + l)C(d + 2), s-0. (2.70) 
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Inserting this into Eq. ( |2.65| ) we immediately recover the zero-temperature result 
(PI)- 

We now seek the leading correction to this. We rewrite I(s) as 

^ = ^Tivi J o i _ e -y/s e J (*-m(y) 



__/ dyyC^e-v/^e-^/'J^D/aCy), (2.71) 

J° n=0 



s (d+5)/2 

where we have employed the geometric series. The Bessel-function integral has an 
elementary form f25fl : 



\j Me -+ = S@£^L ,,72) 
VO + 1 



is the fundamental integral, and the form we want can be written as 

(2.73) 



provided p is a nonnegative integer. (For the application here, this means d is odd, 
but we will be able to analytically continue the final result to arbitrary d.) Then we 
can write I(s) in terms of the series 

(_l)(d+3)/2 °° 

J ( g )= J +5)/2 i ="+ 1 > ( 2 - 74 ) 



where 



^ \ (d+3)/2 



(f) 2 + i-f 



(d-l)/2 



(;) +1 



(2.75) 



We evaluate the sum in Eq. ( p. 74 ) by means of the Euler-Maclaurin summation 
formula, which has the following formal expression: 



5>(Q= / dlF(l) + -[F(oo) + F(l)] 



Here, B n represents the nth Bernoulli number. Since we are considering a ^ 
(the a = term was dealt with in the previous paragraph), the low-temperature 
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limit corresponds to the limit s — > oo. It is easy then to see that F(oo) and all its 
derivatives there vanish. The function F at 1 has the general form 



(d-l)/2 



de / 



e = 1/s -> 0. 



(2.77) 



By examining the various possibilities for odd d, d = 1, d = 3, d = 5, and so on, we 
find the result 



= _(_l)(d-D/2 d M = _ ( _ 1) (d-l)/2 2 (d + l)/2 7r -l/2 r ( « + j 



(2.78) 



Because it is easily seen that before the limit e — > is taken, -F(l) is an even function 
of e, it follows that the odd derivatives of F evaluated at 1 that appear in Eq. (|2.76f) 
vanish. Finally, the integral in Eq. ( |2.76| ) is that of a total derivative: 



dlF(l) 



e \de J 

Thus, the final expression for I(s) is 

1 



Is) 



d 



-F(l). 



[2.79) 



s (d+5)/2 



2 (d-i)/ 27r -i/2 r t + ! \ 



(2.80) 



Note that a choice of analytic continuation to d other than an odd integer has been 
made so as to avoid oscillatory behavior in d. 
We return to Eq. ( |2.65| ). It may be written as 



2Tl d ' 2 



cm 



r(d/2) 2 M 7r d + 1 a d + 1 



a 



(l-d)/2 f 87ra \ 



(d-3)/2 



/ 



v^r(d/2)/ 



a/? 
47ra 



2 -d-l 7r -<i/2-l 



47T 
7 



d+1 



. d 

i i 5 +i 



d+2 ' 



(2.81) 



The correction to the zero-temperature result Eq. fl2.32|) is obtained from 

f T ^° = 

J corr 

= -Tx- dl2 - x Y{dj2 + l)C(d + 2)/r d ~ 2 . 
Thus, the force per unit area in the low-temperature limit has the form 



-^c(2vrn) 

P n=l 



[2.82) 



f 



(d + \)2- d - 2 ix- d l 2 - x 



1 d+2 



T(d/2 + l)((d + 2) 



1 + 



2a 

d+l\~]3 



d+2 



(2.83) 
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of which the d = 1 and d = 2 cases are familiar: 

1 



d= 1 : 
d = 2 : 



./ 



./ 



87ra 3 

r 2 



C(3) 



1 + 4 



/3 3 



7T 



480a 4 



16 a 4 ' 

1 + is 



(2.84a) 
(2.84b) 



These equations are incomplete in that they omit exponentially small terms; for 
example, in the last square bracket, we should add the term 



240 a 



-nP/a 



(2.85) 



Finite temperatures were first discussed by Lifshitz [17]], but then considered more 



fully by Fierz, Sauer, and Mehra |27|, |28|. Hargreaves |29j analyzed the discrepancy 
between the results of Refs. W71 and ^71 . 

2.5 Electromagnetic Casimir Force 

We now turn to the situation originally treated by Casimir: the force between parallel 
conducting plates due to quantum fluctuations in the electromagnetic field. An ele- 
gant procedure, which can be applied to much more complicated geometries, involves 
the introduction of the Green's dyadic, defined as the response function between the 
(classical) electromagnetic field and the polarization source (this formalism is intro- 
duced in Refs. [0,0): 



In the following we will use the Fourier transform of V in frequency: 

r(x, x) = 

which satisfies Maxwell's equations 



^ e -^-*')r(r,r>), 

Z7T 



V x T = 
—V x $ — iuT = iu)15(r — r'). 

The second Green's dyadic appearing here is solenoidal, 

V - $ = 0, 

as is r if a multiple of a 5 function is subtracted: 

v ■ r' = o, r' = r + i5(r - r') 



(2.86) 



(2.87) 



(2.88a) 
(2.88b) 



(2.89) 



(2.90) 
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The system of first-order equations Eqs. ( |2.88a| ), ( |2.88b| ) can be easily converted to 
second-order form: 

(V 2 + u 2 )T' = -V x (V x l)8(r - r'), (2.91a) 
(V 2 + cj 2 )$ = iooV x 15(r-r'). (2.91b) 

The system of equations Eqs. (|2.91a|) , ( |2.91b|) is quite general. We specialize to 
the case of parallel plates by introducing the transverse Fourier transform: 

r'(r, r'; u) = J i^^-^g(z, z'; k, u). (2.92) 

The equations satisfied by the various Cartesian components of g may be easily 
worked out once it is recognized that 

[V x (V x l)} tJ = V. 4 V, - <%V 2 . (2.93) 

In terms of the Fourier transforms, these equations are 

f |1 - k 2 + u 2 ) g zz = -k 2 5(z - z'), (2.94a) 



k 2 + lu 2 ) g zx = -ik x ^8{z - z>), (2.94b) 



oz l } OZ 



k 2 + LU 2 jg zy = -ik y —5(z-z'), (2.94c) 
°~ k 2 + lu 2 ] g xx = i-k 2 + 5{z - z'), (2.94d) 



\dz 2 



OZ z J \ OZ 

[^~k 2 + ^ g yy = (-k 2 x + ^ 8(z - z% (2.94e) 

(J? ~ k2 + 9xy = kxky5{ - Z ~ Zl) ' (2 ' 94f) 

We solve these equations subject to the boundary condition that the transverse 
components of the electric field vanish on the conducting surfaces, that is, 

n x r'| 2=0 , a = 0, (2.95) 

where n is the normal to the surface. That means any x or y components vanish at 
z = or at z = a. Therefore, g xy is particularly simple. By the standard discontinuity 
method, we immediately find [cf. Eq. ( [2.191) 1 

k k 

9xy = g yx = . J \ (ss), (2.96) 

A Sill ACL 



19 



where 

(ss) = sin Xz < sin X(z > — a). (2.97) 
To find g xx we simply subtract a 5 function: 

g' xx = g xx -5(z-z'). (2.98) 

Then, we again find at once 

A sin Aa 

and similarly 

I'** = S < SS >- (2 ' 100) 

To determine the boundary condition on g zz , we recall the solenoidal condition on T', 
Eq. (|2.9U|), which implies that 



d 

—9zz 



dz 



= 0. (2.101) 

,z=0, a 



This then leads straightforwardly to the conclusion 

^ = -Y-^-( cc )' (2.102) 
A sin Aa 

where 

(cc) = cos Xz < cos A(z> — a). (2.103) 

The remaining components have the property that the functions are discontinuous, 
while, apart from a 5 function, their derivatives are continuous: 

9zx = —^-(c8), (2.104) 
sin Aa 

9zy = -^-(cs), (2.105) 
sin Aa 

where, because of the analogue of Eq. ( |2.101| ), or because the components must vanish 
when z' = 0, a, 

| cos Xz sin \(z> -a) z<z> 
I sin Xz cos X{z — a), z > z . 



20 



Similarly, 



sin Xa ' 



sc 



9yz 



^>ky 

sin Xa 



sc 



where 



[sc) 



sin Xz cos X(z' — a), z < z', 
cos Az' sin A(z — a), z > z'. 



(2.107) 



(2.108) 



(2.109) 



The relation between Eqs. flUMp , ( gTlOg ) and Eqs. (glQ?|) , (|2~T08|) just reflects the 
symmetry 



r'(r,r') = r' r (r',r), 
where T denotes transposition, which means that 

g lj (z,z';'k) = gji(z',z; -k). 
The normal-normal component of the electromagnetic stress tensor is 

T ZZ = 1 -(E 2 + H 2 )-(E 2 Z +H 2 ). 

The vacuum expectation value is obtained by the replacements 

(E(z)E(a/)) = -T(x,x'), 



(2.110) 
(2.111) 

(2.112) 
(2.113a) 



(H(x)H(x')) 



1 1 



% or 



v x r(x,x) x V'. 



In terms of the Fourier transforms, we have 



[-(oo 2 - k 2 )g zz + (uj 2 - k 2 y )g xx + (oo 2 - k\ 



+ *ky y dz 9yz dz ,9zy j + ik x ^g^J 
~\~ k x ky{cj X y -\- g-yx) ~\~ ~~q^ {dxx 4" 9yy) 



(2.113b) 



(2.114) 



When the appropriate Green's functions are inserted into the above, enormous sim- 
plification occurs on the surface, and we are left with 



{T zz )\ 



iX cot Xa, 



I z=z'=0,a 

which indeed is twice the scalar result ( [2.21| ), as claimed at the end of Sec. pTT 



2.115) 
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2.6 Fermionic Casimir Force 



We conclude this section with a discussion of the force on parallel surfaces due to 
fluctuations in a massless Dirac field. For this simple geometry, the primary distinc- 
tion between this case and what has gone before lies in the boundary conditions. The 
boundary conditions appropriate to the Dirac equation are the so-called bag-model 
boundary conditions. That is, if n M represents an outward normal at a boundary 
surface, the condition on the Dirac field ip there is 

(1 + m- 7)^ = 0. (2.116) 

For the situation of parallel plates at z = and z = a, this means 

(l T i 7 3 )^ = (2.117) 

at z = and z = a, respectively. In the following, we will choose a representation of 
the Dirac matrices in which 275 is diagonal, in 2 x 2 block form, 




(2.118) 



while 

from which the explicit form of all the other Dirac matrices follow from 7 = i^^^a . 



2.6.1 Summing modes 



It is easiest, but not rigorous, to sum modes as in Sec. |2TT| . We introduce a Fourier 
transform in time and the transverse spatial directions, 

duJ „-farf f W „ik.x„ 



^) = J ^ e J J^ e 1 ^;k,u;), (2.120) 

so that the Dirac equation for a massless fermion —i^dip = becomes (if k lies along 
the negative y axis) 

d \ 

-Lu^i—)u±±kv± = 0, (2.121a) 
oz J 

±ku±+ f-u±i-^j v± = 0. (2.121b) 

where the subscripts indicate the eigenvalues of 275 and u and v are eigenvectors of 
a 3 with eigenvalue +1 or —1, respectively. This system of equations is to be solved 
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to the boundary conditions ( |2.117| ), or 



u + + U-\ z=0 = 0, 


(2.122a) 


v+ - U-|z=0 = 0, 


(2.122b) 


u + - U-\ z=a = 0, 


(2.122c) 


v + + v^\ z=a = 0. 


(2.122d) 



The solution is straightforward. Each component satisfies 

^ = 0, 



where A 2 



to 



— A< 

dz 2 



k 2 , so that the components are expressed as follows: 



U + + M_ 
V + — t>_ 
U + — M_ 

V+ + u_ 



A sin Xz, 
B sin Xz, 
C sin A(z — a), 
D sin A (z — a). 



(2.123) 



(2.124a) 
(2.124b) 
(2.124c) 
(2.124d) 



Inserting these into the Dirac equation (|2.121a| ) and (|2.121b|) , we find, first, an eigen- 
value condition on A: 



or 



cos Xa = 0, 



1 



(2.125) 



\a = ( n +^)7T, (2.126) 

where n is an integer. We then find two independent solutions for the coefficients: 

(2.127a) 
(2.127b) 

(2.127c) 



A^0, 

B = 0, 

iuj 



C 



and 



A = 0, 



-l) n A, 
-ITA, 



C 
D 



iX 

to 

~x { 



-l) n B, 
-l) n B. 



(2.127d) 

(2.128a) 
(2.128b) 

(2.128c) 
(2.128d) 
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Thus, when we compute the zero-point energy, we must sum over odd integers, noting 
that there are two modes, and remembering the characteristic minus sign for fermions: 
Instead of Eq. (|2.2|) , the Casimir energy is 



2 ^ (2vr) 2 V a 2 

1 1 f°° dt . _ 3/2 t (n+l/2) 2 w 2 /a 2 



2JnA7T ^ In t 
v n=0 JU 



-r (-?) f 



(ra+ 1/2) 3 tt 3 

x ' n=0 

= -^g«-3), (2.129) 

which is | x 2 times the scalar result (|2.8[) because C( — 3) = — S4/4 = 1/120. (The 
factor of 2 refers to the two spin modes of the fermion.) 

2.6.2 Green's Function Method 

Again, a more controlled calculation starts from the equation satisfied by the Dirac 
Green's function, 

-y-dG{x,x') = S(x-x'), (2.130) 

% 

subject to the boundary condition 

(l+m- 7 )G 

We introduce a reduced, Fourier-transformed, Green's function 



= 0. (2.131) 

2=0, a 



which satisfies 



G{x,x>) = j j ^ 2 ^-^g{z,z'-M, (2.132) 

f + 7 . k + g{Zj z>) = s{z _ z ,y (2 133) 

Introducing the representation for the gamma matrices given above, we find that the 
components of g corresponding to the +1 or —1 eigenvalues of 175, 
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satisfy the coupled set of equations 



, d , 

-u ± a ■ k T *o" ^ ) fci = 
d 

-u ± cr • k =F ) 9±t = Ti5(z 



We then resolve each of these components into eigenvectors of cr 3 : 



9±± 



„(+) „(-) 
w (+) ,,(-) 



and similarly for g± T . These components satisfy the coupled equations 

d 



-oj =F i— u\l ± fcf 
<7.; 



: : 



±Jfcu£2 + 







,(±) 



<9 



(+) 



(-) 

±=F 



, (-) 



0. 

o. 

0. 
0. 

TiS(z - z'), 



(2.135a) 
(2.135b) 



(2.136) 



which aside from the inhomogeneous terms are replicas of ( 2.121a]) and 
These equations are to be solved subject to the boundary conditions 

u {±) 
"++ 

(±) 



(2.137a) 
(2.137b) 
(2.137c) 
(2.137d) 
(2.137e) 
(2.137f) 
( |2.121b|) . 



u 



a 



u 



(±) 

+- 

(±) 

+- 

; (±) 



-u^l\ z=a = 0, 


(2.138a) 


+ u^l\ z=0 = 0, 


(2.138b) 


-u {± l\ z=a = 0, 


(2.138c) 


+ ^|, =0 = 0, 


(2.138d) 


+ w^U=« = o, 


(2.138e) 


-v^\ z=0 = 0, 


(2.138f) 


+ ^ ± _ ) |* =a = 0, 


(2.138g) 


-^1^0 = 0. 


(2.138h) 
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Again, the solution is straightforward. We find 



(+)* 



(-)* 
v + _f = v 



(-) 



2 cos Aa 



cos \(z + z' — a) + ^ sin A(z + z' — a) 
A 



v (+) _ 



,(+) = ..H 
ifc 1 



2A cos Aa 



sin A(z + z' — a), 



(2.139a) 



(2.139b) 



(+) 



(+)* 



.(-)* 



,(-) 



2 cos Aa 



[e(z — z') cos \{z > — z< — a) 



— sin A(^> — z < — a) 
A 



,(+) 



v : 
ik 



(+) 



2A cos Aa 



sin X(z > — z< — a), 



where 



1 if z > z', 
-1 if z < z'. 



(2.139c) 



(2.139d) 



(2.140) 



We now insert these Green's functions into the vacuum expectation value of the 
energy-momentum tensor. The latter is 



(2.141) 



We take the vacuum expectation value by the replacement 

^7° - -G, 

% 



(2.142) 



where G is the fermionic Green's function computed above. Because we are interested 
in the limit as x' — > x we can ignore the Lagrangian term in the energy-momentum 
tensor, leaving us with in the transform space 



(T 



Tr^d^z') 
tr a 3 (g- + + 



1 
2 

29z 



5^+«S2 -(»3+«S>] 



(2.143) 
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When we insert the solution found above, Eqs. ( p.l39a| )-( |2.139d| ), we obtain 

( T 33) =2 ,|_ ReM (+) 
OZ 



(2.144) 
(2.145) 



Carrying out the differentiation and setting z = z' we find 

(T 33 ) = zAtanAa, 

where again we ignore the 5-function term. [Compare the scalar result ( 2.21|) .] 
We now follow the same procedure given in Sec. [2.2| : The force per unit area is 

d 2 k 



f 



(2vr) 2 
d 2 k 



duo 
2tt' 



-iX tan Xa 



— ptanh pa 
2% 



2tt 2 



dpp 3 



1 



e 2pa + I 



(2.146) 



As in Eq. (|2.23 ) we omit the 1 in the last square bracket: The same term is present 
in the vacuum energy outside the plates, so cancels out when we compute the discon- 
tinuity across the plates. We are left with, then, 

"°° dxx 3 



f 



1 



But 



so here we find 



16n 2 a 4 
00 X s " 1 dx 



e x + 1 



o 



e x + 1 



ft 



1 - 



7tt 2 



1920a 4 ' 



(2.147) 



(2.148) 



(2.149) 



which is, indeed, % times the scalar force given in Eq. (|2.9|). 



The effect of fermion fluctuations was first investigated by Ken Johnson [32 



in 



connection with the bag model |33 



3 Casimir Effect in Dielectrics 



The formalism given in Sec. |2.5| can be readily extended to dielectric bodies |30[. The 
starting point is the effective action in the presence of an external polarization source 
P: 

W = J (dx)[P ■ (-A - V0) + eE ■ (-A - V0) - H • (V x A) + ^(H 2 - eE 2 )], 

(3.1) 
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which, upon variation with respect to H, E, A, and 0, yields the appropriate 
Maxwell's equations. Thus, because W is stationary with respect to these field vari- 
ations, its response to a change in dielectric constant is explicit: 



S e W = J (dx)5e^E 2 . 



(3.2) 



Comparison of iS e W with the second iteration of the source term in the vacuum 
persistence amplitude, 



jw 



1 

+ 2 



i J (dx)E ■ P + 



allows us to identify the effective product of polarization sources, 

iP(x)P(x')| eff = I5e5(x-x'). 
Thus, the numerical value of the action according to Eq. (|2.86|) , 

W=^J (dx)P(x) ■ E{x) = ^J (dx)(dx')P(x) -T{x,x') • P(x') 



(3.3) 



(3.4) 



(3.5) 



implies the following change in the action when the dielectric constant is varied 
slightly, 



5W = -- 



(dx)Se(x)T kk (x,x), 



(3.6) 



which, in view of Eq. (|3.2| ), is equivalent to the vacuum-expectation- value replacement 

(gmp . 

For the geometry of parallel dielectric slabs, like that considered in Sec. 0, but 
with dielectric materials in the three regions, 



ei, z < 0, 
r( : s = < e 3 , < z < a, 
e 2 , a < z, 



(3.7) 



the components of the Green's dyadics may be expressed in terms of the TE (trans- 
verse electric or H) modes and the TM (transverse magnetic or E) modes, given by 
the scalar Green's functions satisfying 



d 2 



dz 2 



+ k 2 -u 2 t )g E {z,z') = 6{z-z'), 



(3.8a) 
(3.8b) 
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where, quite generally, e = e(z), e' = e(z'). The nonzero components of the Fourier 
transform g given by Eq. ( |2.92|) are easily found to be (in the coordinate system where 
k lies along the +x axis) 







1 d 1 d H 
+ e~d~z7'lh' 9 ' 


(3.9a) 


9xx 


= —5(z - z) 


n 

yyy 


- u 2 a E 




C3 9bl 


9zz 


= --5(z - z') 




(3.9c) 




ee' 


9xz 


. k d H 
ee' oz 




(3.9d) 


9zx 


. k d H 
ee' oz' 




(3.9e) 



The trace required in the change of the action ( |3.6| ) is obtained by taking the limit 
z' — > z, and consequently omitting delta functions: 



9kk = I oj 2 g E + —g H + - — -— g H 
• ee' e oz e' oz' 



(3.10) 



This appears in the change of the energy when the second interface is displaced by 
an amount Sa, 

Se(z) = — <5a(e 2 — e 3 )5(z — a), (3-H) 
namely (A is the transverse area) 



/ ^^ dz6 < z )9kk(z,z';k,u). (3.12) 



SE i f du (dk) 
~A = 2 



Because g E , g H and \§i]j^g H are all continuous, while ee' is not, we interpret the 
trace of g in Eq. ( |3.12|) symmetrically; we let z and z' approach the interface from 
opposite sides, so the term ^pg H — > -^g H - Subsequently, we may evaluate the 
Green's function on a single side of the interface. In terms of the notation 

n 2 = k 2 -uj 2 e, (3.13) 

which is positive after a complex frequency rotation is performed (it is automatically 
positive for finite temperature), the electric Green's function is in the region z, z' > a, 

g E ( Z , Z') = — (e- K ^ Z - Z '\ + re -^{z+z'-2a)\ ^ 



2k 2 

where the reflection coefficient is 

K 2 - K 3 , 4k 2 K 3 



d~\ (3.15) 
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with the denominator here being 

d= K3 + Kl *»±*V "»«-1. (3.16) 

«3 - «1 «3 - ^2 

The magnetic Green's function g ff has the same form but with the replacement 

k — ► n/e = (3-17) 

except in the exponentials; the corresponding denominator is denoted by d'. [It is 
easy to see that g E reduces to Eq. Q2.24 ) when r = — 1; the results in Sec. 2.5 follow 
from Eqs. ( |3.9a|) -( j3.9e|) in the coordinate system adopted here.] 

Evaluating these Green's functions just outside the interface, we find for the force 
on the surface per unit 

f = \J |^^r { [«a - ^2 + 2^d- x ] + [k 3 - k 2 + 2^d'~ l ] } , (3.18) 

where the first bracket comes from the E modes, and the second from the H modes. 
The first term in each bracket, which does not make reference to the separation a 
between the surfaces, is seen to be a change in the volume energy of the system. 
These terms correspond to the electromagnetic energy required to replace medium 2 
by medium 3 in the displacement volume. They constitute the so-called bulk energy 
contribution. The remaining terms are the Casimir force. We rewrite the latter by 
making a complex rotation in the frequency, 

w->zC, Creal. (3.19) 

This gives for the force per unit area at zero temperature 

-I /*oo /*oo 

fZ^r = -7- 2 d( dk 2 2k 3 (cT 1 + d'- 1 ) , K 2 = fc 2 + eC 2 . (3.20) 

Jo JO 

From this, we can obtain the finite temperature expression immediately by the sub- 
stitution 

C 2 C„ 2 = taV/,3 2 , (3.21) 



n=0 



the prime being a reminder to count the n = term with half weight. These results 
agree with those of Lifshitz et al. |T%|| . 

Note that the same result ( 3.20|) may be easily rederived by computing the normal- 



normal component of the stress tensor on the surface, T zz , provided two constant 
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stresses are removed, terms which would be present if either medium filled all space. 
The difference between these two constant stresses, 



T. 



vol 



(dk±) duj , 
(2tt) 2 2^ 



«2 - «3 



(3.23) 



precisely corresponds to the deleted volume energy in the previous calculation. 

Various applications can be made from this general formula (|3.20|) . In particular, 
if we set the intermediate material to be vacuum, 63 = 1, and set K\ — K2 — 00, 
/"vi — Kj<2 = 0, we recover the Casimir force between parallel, perfectly conducting 



plates, Eq. ( |2.10| ), including the correct temperature dependence, given by Eq. (|2.51| ) 
with d = 2 and an additional factor of two representing the two polarizations. An- 
other interesting result, important for the recent experiments |]20|| , is the correction 
for an imperfect conductor, where for frequencies above the infrared, an adequate 
representation for the dielectric constant is [|4[ 



e(u;) 



where the plasma frequency is 



to' 



1 



4vre 2 iV 



m 



(3.24) 



(3.25) 



where e and m are the charge and mass of the electron, and N is the number density 
of free electrons in the conductor. A simple calculation shows, at zero temperature 



/ 



7T 



240a 4 



1 



3t/tt ea 



m\V 2 ' 

n) 



(3.26) 



Now suppose the central slab consists of a tenuous medium and the surrounding 
medium is vacuum, so that the dielectric constant in the slab differs only slightly 
from unity, 



e-K 1. 



(3.27) 



Then, with a simple change of variable, 



K 



(P, 



(3.28) 



4 Although we have used rationalized Heaviside-Lorentz units in our electromagnetic action for- 
malism, that is without effect, in that the one-loop Casimir effect is independent of electromagnetic 
units. For considerations where the electric charge and polarizability appear, it seems more conve- 
nient to use unrationalized Gaussian units. 
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we can recast the Lifshitz formula ( |3.20| ) into the form 

fVc 3 r^KC)-l] 2 [(2p 2 -l) 2 + l]e- 2 ^. (3.29) 
Jo Ji V 

If the separation of the surfaces is large compared to the characteristic wavelength 
characterizing e, a( c ^> 1, we can disregard the frequency dependence of the dielectric 
constant, and we find 

23(e-l) 2 , . 

f - -ekd- (3 - 30) 

For short distances, a( c <C 1, the approximation is 

f^-^-J ^C(e(C) " I) 2 - (3.31) 



32vr 2 a 3 .„ 

These formulas are identical with the well-known forces found for the complementary 
geometry in [pOj]. 

Now we wish to obtain these results from the sum of van der Waals forces, derivable 
from a potential of the form 

V = - — . (3.32) 
We do this by computing the energy (jV = density of molecules) 

E = -\ B #[iz[*/j { ^ )(iA) w -^L__ . (3.33) 

If we disregard the infinite self-interaction terms (analogous to dropping the volume 
energy terms in the Casimir calculation), we get 

f- dE - 2 * BN2 1 (3 34) 

J da A (2 -7) (3 -7) a?" 3 ' 1 " ; 

So then, upon comparison with ( |3.30| ), we set 7 = 7 and in terms of the polarizability, 

a = i^' (3 - 35) 



we find 



23 

B = —a 2 , (3.36) 

47T 

or, equivalently, we recover the retarded dispersion potential [[|, 

23 a 2 , 
V = (3.37) 
47r r' 

whereas for short distances we recover from Eq. ( |3.31| ) the London potential 0, 

3 1 f°° 

V = —- rfC«(C) 2 , (3-38) 
which are the quantititive forms of Eqs. (|1.3|) and ( |1.2| ), respectively. 
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4 Casimir Effect on a Sphere 



4.1 Perfectly Conducting Spherical Shell 

The zero-point fluctuations due to parallel plates, either conducting or insulating, 
give rise to an attractive force, which seems intuitively understandable in view of 
the close connection with the attractive van der Waals interactions. However, one's 
intuition fails when more complicated geometries are considered. 

In 1956 Casimir proposed that the zero-point force could be the Poincare stress 



stabilizing a semiclassical model of an electron |35| . For definiteness, take a naive 
model of an electron as a perfectly conducting shell of radius a carrying a total 
charge e. The Coulomb repulsion must be balanced by some attractive force; Casimir 
proposed that that could be the vacuum fluctuation energy, so that the effective 
energy of the configuration would be 

e 2 Z 

E= e - he, (4.1) 

2a a 

where the Casimir energy is characterized by a pure number Z. The would open the 
way for a semiclassical calculation of the fine-structure constant, for stability results 
if E = or 

e 2 

a = — = 2Z. (4.2) 

he 



Unfortunately as Tim Boyer was to discover a decade later ||36[| , the Casimir force 
in this case is repulsive, Z = —0.04618. The sign results from delicate cancellations 
between interior and exterior modes, and between TE and TM modes, so it appears 
impossible to predict the sign a priori. 

Boyer's calculation was rather complicated, involving finding the zeroes of Bessel 
functions. In 1978 we found a simpler approach, based on Green's functions [3~I], 



which I will describe here. In particular the Green's dyadic formalism of Sec. 2.5 may 



be used, except now the modes must be described by vector spherical harmonics, 



defined by [34, 37 



X lm =[l(l + l)}^ 2 LY lm (0,<P), (4.3) 
where L is the orbital angular momentum operator, 



L = -r x V. (4.4) 

i 



The vector spherical harmonics satify the orthonormality condition 

J dOX*/ m / • X; m = 5ui5 mm i, (4.5) 



33 



as well as the sum rule 

1 

" \X,„(9.(h)\ 2 = - 

47T 



£ \*ueA)\ 2 = 2± ^- (4.6) 



m=—l 



The divergenceless dyadics V' and <fr may be expanded in terms of vector spherical 
harmonics as 

r ' = ( fiXim + ^ v x 9iXim ) ' ( 4 - 7a ) 

$ = ^ (&X Zm - -V x /,X Jm ) , (4.7b) 

where the second suppressed tensor index is carried by the coefficient functions fi, gi, 

fh 9i- 

Inserting this expansion into the first-order equations ( j2.88a| ) and (|2.88b[ ), and 
using the properties of the vector spherical harmonics, we straightforwardly find 
that the Green's dyadic may be expressed in terms of two scalar Green's functions, 
the electric and the magnetic: 



r(r,r';u,) = X> 2 F,(r, r')X lm ({l)X* m (a) 



Im 

+ 5-function terms, (4.8) 

where the expression "5-function terms" refers to terms proportional to spatial delta 
functions. These terms may be omitted, as we are interested in the limit in which 
the two spatial points approach coincidence. These scalar Green's functions satisfy 
the differential equation 

l -^r - l «±H + A I F ^ r ') \ = -U(r - r>) (4 9) 

rdr* r2 + J {G^r') j r 2 h K ^> 

subject to the boundary conditions that they be finite at the origin (the center of 
the sphere), which picks out ji there, and correspond to outgoing spherical waves at 
infinity, which picks out h\ \ . On the surface of the sphere, we must have 



^(0,0 = 0, ^rG^r') 



0. (4.10) 



r=a 



The result is that 

Fij~~ l {Ft 



Cl] G° + \ G J\, (4.11) 
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where G° is the vacuum Green's function, 

G?(r,r')=ikj l (kr < )h\ 1 \kr > ), (4.12) 
and in the interior and the exterior of the sphere respectively, 

r,r' < a : j G J } = -A GtF ikji (kr)ji (kr 1 ) , (4.13a) 

r,r' > a : j ^ | = -B G}F ikhf ] \kr)hf' '(At'), (4.13b) 

where the coefficients are 

^-S* (414b) 

From the electromagnetic energy density we may derive the following formula for the 
energy of the system 

E = |(dr)i /^e-^')^{A; 2 [FKr,rO + GKr,rO]X /m (fi).XL(fi') 



- V x X lm (Q) ■ [F,(r,r') + G,(r,r')] . X^^) x V} 



(4.15) 



Note here that the vacuum term in the Green's functions has been removed, since 
that corresponds to the zero-point energy that would be present in this formalism 
if no bounding surface were present. Here we are putting the two spatial points 
coincident, while we leave a temporal separation, t — t — t', which is only to be set 
equal to zero at the end of the calculation, and therefore serves as a kind of regulator. 
The integration over the solid angle and the sum on m may be easily carried out, 
and the radial integral over Bessel functions is simply done using recurrence relations. 
After performing a complex frequency rotation, we obtain the following compact form 
for the Casimir energy: 

1 °° 1 f°° d 

E = -^T a B 2i + !>2 J_J ye "" X di ln(1 " A?) ' ( " 6) 

where 

A, = [s^dix)}' (4.17) 
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is written in terms of Ricatti-Bessel functions of imaginary argument, 



si(x) 



ei(x) 



TlX 




h+i/2{x), 




2x 

—K, 
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1+1/2 



X) 



Here, as a result of the Euclidean rotation, 



1 lr 

x = \y \, y = -ka is real, as is e = > 0. 

i i a 



(4.18) 



(4.19) 



The same formula may be derived by computing the stress on the surface through 



use of the stress tensor 31 



A very rapidly convergent evaluation of this formula can be achieved by using the 
uniform asymptotic expansions for the Bessel functions: 



s t (x) 



edx 



4/2 



2(l + z 2 ) 1/2 
z 1 ' 2 _ 

( 1 + z 2)l/2 e 



VI] 



u k (t) 



k=l 
oo 



l) k u k (t) 



k=l 



, / — > OO, 



where 



x = uz,i/ = l + 1/2, t=(l + z 2 y 1/2 , rj = (1 + z 2 ) 1/2 + In 



1 + (1 + Z 2 )V2 



(4.20a) 
(4.20b) 

(4.21) 



and the Ufe(t) are polynomials in t of definite parity and of order 3k f38 |. If we now 
write 



oo 



we easily find 



J(Z,0)~— , /^oo. 



(4.22) 



(4.23) 



In order to obtain a finite sum, therefore, we must keep e ^ until the end of the 
calculation. By adding and subtracting the leading approximation to the logarithm, 
we can write 



J(/,e)=i2, + 5 { (e), 



where 



Ri 



2tt 



dz 



(2/ + l)^ln(l-Af) +(1^ 



J(Z,0)- 



32' 



(4.24) 



(4.25) 
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and 



Site) 



1 

47T 



z dz e % 



d 



dz{l + z 2 Y 

By use of the Euler-Maclaurin sum formula Q2.76 ), we can work out the sum 



1=1 



3 

3~2 : 



(4.26) 



(4.27) 



precisely the negative of the value of a single term at e = 0! The sum of the remainder, 
Ri, is easily evaluated numerically, and changes this result by less than 2%. Thus 
the result for the Casimir energy for a spherical conducting shell is found to be 



E 



0.092353 
2a ' 



(4.28) 



This agrees with the result found in 1968 by Boyer |pq|, evaluated more precisely by 



Davies [[$9j , and confirmed by a completely different method by Balian and Duplantier 
in 1978 1 40| . Recently, this result has been reconfirmed, using a zeta function method, 
by Leseduarte and Romeo [|T . A reconsideration using direct mode summation has 
appeared very recently [02 . 



4.2 Fermion Fluctuations 



The corresponding calculation for a massless spin-1/2 particle subject to bag model 
boundary conditions ( |2.116|) on a spherical surface, 



(l + m- 7 )G 



0. 



(4.29) 



was carried out by Ken Johnson |43[ and by Milton |44[ . The result is also a repulsive 
stress, of less than one-half the magnitude of the electromagnetic result. (Recall that 
for parallel plates, the reduction factor was 7/8.) 

In this case we wish to solve the Green's function equation 



7-<9 ) G(x, x') = 5(x — x 1 ) 



(4.30) 



subject to the boundary condition ( |4.29|) . In the same representation for the gamma 
matrices used before, this may be easily achieved in terms of the total angular mo- 
mentum eigenstates ( J = L + (l/2)cr): 



7 l=J±l/2 
' JM 



(n) 



1 + 1/2TM 



21 + 1 



1/2 



1 + 1/2 ±M 
21 + 1 



1/2 



Y, 



IM+1/2 



(n)|->. 



(4.31) 
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These satisfy 

a ■ rZ^f 1/2 = Z 1 ^ 112 - (4.32) 

Once the Green's function is found, it can be used in the usual way to compute the 
vacuum expectation value of the stress tensor, which in the Dirac case is given by 

= -^ 7 °- (W + Y-dA V> + sTC, (4.33) 



2 ' 2 V i i 
£ being the fermionic Lagrange function, which leads directly to 

1 d 

T rr = oTr tr 7 ■ rG(x,x) 

2 or 



(4.34) 



The discontinuity of the stress tensor across the surface of the sphere gives the energy 
according to 

g 

4na 2 [T rr (a-) - T rr (a+)] = --E(a). (4.35) 



The result of a quite straightforward calculation (the details are given in Ref. |45| ) 
gives the result for the sum of exterior and exterior modes, again, in terms of modified 
spherical Bessel functions: 

2 °° f°° d 

E = YVZ + l) / dxxcosxe — ln^ + e 2 l+1 ){s 2 l + s 2 l+1 )\. (4.36) 

7ia Jq ax 

This may again be numerically evaluated through use of the uniform asymptotic 
approximants, with the result 

E = °J™. (4.37) 
a 

Somewhat less precison was obtained because, in this case, the leading uniform asymp- 
totic approximation vanished. 



5 Infinite Circular Cylinder 

Since parallel plates yield an attractive Casimir force, and a sphere has a repulsive 
stress, one might guess that for a cylinder a zero stress results. The situation is not so 



simple. The first calculation was carried out in 1981 by DeRaad and Milton ||46|| . The 
electro dynamic result turns out to be attractive but with rather small magnitude. 

Consider a right circular perfectly conducting cylinder of infinite length and radius 
a. We compute the Casimir energy using the above Green's dyadic formalism, adapted 



38 



to this cylindrical basis. The necessary information about vector spherical harmonics 
in this case is given in Stratton p7f. The results for the Green's dyadics are 



r(r,r') 



E x f dk 
I 2vr 



m=—oo 



-— MM'*(d m -k 2 )F m (r } r') 



H ^NN'*G m (r, r') 



10 

oo 



X m k( , Z )Xmk(°', Z '), 



<fr(r, r') 



Ef dk 
I 2vr 



VTL— — OO 



iMN'*G m (r, r') 



- -NM"F m (r, r') 



where the nonradial eigenfunctions are 

Xmk(0,z) : 



Xmk(°, Z )X*mk(°', Z ' 



^imO ^ikz 



the vector differential operators M and N are 



^im ~d ^ T d A mk 

M = f 0— , N = hk- zd m , 

r or or r 

and the scalar differential operator d m is 

1 d d m 2 

The scalar Green's functions in the interior and exterior regions are 

H'JXa) 



r,r < a : 



-G m (r,r') 



,/ral r i r ) — 



J m (Ar<)# m (Ar>) 



-J m (\r)J m (\r') 



171 



+ ^&,r>), 



2m \r 



± 



J m (Xa) 

Y j~|w|~,/|m| 

2\m\ a 2 H ' 



(5.1a) 



(5.1b) 



(5.2) 



(5.3) 



(5.4) 



(5.5a) 

(5.5b) 
(5.5c) 
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r, r > a : 
a; 2 m[ ' j 2A 2 



J m (Ar<)# m (Ar>) 



1 



a; 



-G m (r,r') 



27T 

2A 2 



J m (Ar<)iJ m (Ar>) - 



J m (Aa) 
H m (Xa) 



H m {\r)H m {\r') 



H m (Xr)H m (Xr') 



r, r') 



2 TO 



1 fl2|m| TO^ 

r> / 2|m| rHr'H' m ' 



(5.6a) 

(5.6b) 
(5.6c) 



where H m = Hm is the Hankel function of the first kind, and A 2 = uo 2 — k 2 . Although 

G F 

Q ' are not determined, they do not contribute to physical quantities. 
The result for the force per unit area is (z = \a) 



f = T rr (a_) - T rr (a + ) 



.(57) 



Here, because convergence is more subtle than in the previous cases, we have inserted 
a frequency cutoff function ip{uj) that has the properties 

^(0) = 1, 

ip(uj) ~ - — — , la; I — > oo, 



\u> 



ip{oj) real for cu real or imaginary. 
A way of satisfying these conditions is to take 



m = E 



+ 



a 2 oj 2 — jif ' a 2 oj 2 — m 



(5.8) 



(5.9) 



with the conditions on the residues and poles: 

^ReOj = 0, 

i 



(5.10a) 
(5.10b) 



Of course, the poles are to recede to infinity. When we rotate the contour to imaginary 
frequencies, the contribution of these poles is necessary to achieve a real result. 
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Once again we use the uniform asymptotic expansion to extract the leading be- 
havior. We evaluate the m sum for the leading term by using the identity (|2.42|) , 
or 



oo 1 
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m=l 



m^p — 



V~ z2 p 



1 — cotli7n 



TT 



The result for the energy per unit length is 



£ = %a 2 f 



V~ z2 p 



iS + R + Ro) 



(5.11) 



(5.12) 



where 



1 f 00 

S — - I dr 7r(coth7rr — 1) 

1 / d 7 d 2 1 d 3 \ f°° , [ ix / nr \ 1 
+ - 5— H — r H — / dr — coth 1 

2 V dp 2dp* 2dp3) J [^p V Jp J r 



1, n 5 
— In2vre + -, 

2 8' 



i?0 



°° POO ( 

-4 ^ / r In [l - {r{I m {r)K m {r))'f] + 

m=l ^ ^ 



-0.0437 



4(m 2 + r 2 ) 3 

2^°° o rdr jln [l - {r{h{r)K {r))') 2 ] + ^} - ~ 



- In e + 0.6785. 



(5.13a) 



(5.13b) 



(5.13c) 



Adding these numbers, we see that the infrared singularity e cancels, and we obtain 
an attractive result, 



£ 



0.01356 



(5.14) 



Very recently, this result has been confirmed by two independent calculations. 
First, Gosdzinsky and Romeo obtained the same answer, to eight significant figures, 
using a zeta-function technique ||47|| . Shortly thereafter, an earlier calculation by 
Nesterenko was corrected to yield the same answer | 4"8]| . This later method is based 
not on the Green's function, but rather an analytic evaluation of the sum of zero-point 
energies through the formal formula 



(5.15) 
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where f n is a function, the zeroes of which are the mode frequencies of the system, 
and where C is a contour initially chosen to encircle those zeroes. The divergences 
there are regulated by a zeta function technique; the results, both analytically and 
numerical, are identical to those reported here and derived in Ref. |16 



6 Casimir Effect on a D- dimensional Sphere 

Because of the rather mysterious dependence of the sign and magnitude of the Casimir 
stress on the topology and dimensionality of the bounding geometry, we have carried 
out calculation of TE and TM modes bounded by a spherical shell in D spatial 
dimensions. We first consider massless scalar modes satisfying Dirichlet boundary 
conditions on the surface, which are equivalent to electromagnetic TE modes. Again 
we calculate the vacuum expectation value of the stress on the surface from the 
Green's function. 

The Green's function G(x, t; x', t') satisfies the inhomogeneous Klein-Gordon equa- 
tion (gig ), or 



~ V ') G(X ' t; ^ = 6(D){ * ~ X ' )5(t " *' } ' (6 - 1} 

where V 2 is the Laplacian in D dimensions. We will solve the above Green's function 
equation by dividing space into two regions, region I, the interior of a sphere of radius 
a and region II, the exterior of the sphere. On the sphere we will impose Dirichlet 
boundary conditions 

G(x,t;x , ,t') | |x|=a =0. (6.2) 

In addition, in region I we will require that G be finite at the origin x = and in region 
II we will require that G satisfy out going- wave boundary conditions at |x| = oo. 

The radial Casimir force per unit area / on the sphere is obtained from the radial- 
radial component of the vacuum expectation value of the stress-energy tensor ||31|| : 



/=(0|2j:-2^|0> | r=o . (6.3) 

To calculate / we exploit the connection between the vacuum expectation value of the 
stress-energy tensor T M1/ (x, t) and the Green's function at equal times G(x, £; x', t): 



"5 



d d d d 

tttt-;G(x, t; x , t)- m - — — G(x, t; x , t) out 
or or or or 



(6.4) 



x=x', x =a 



To evaluate the expression in (|6.4j) it is necessary to solve the Green's function 
equation (|6.1|) . We begin by taking the time Fourier transform of G: 

/CO 
dte Mt-t') G ^ t]X ' jt 'y ( 6 .5) 
-co 
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The differential equation satisfied by is 

-(^ 2 + V 2 )^(x;x') = ^(x-x'). 



(6.6) 



To solve this equation we introduce polar coordinates and seek a solution that 
has cylindrical symmetry; i.e., we seek a solution that is a function only of the two 
variables r = |x| and 8, the angle between x and x' so that x • x' = rr' cos 9. In terms 
of these polar variables ( |6.6| ) becomes 

d 2 



dr 2 



D — 1 d sin 2 ~ D 6 8 D 2 „ d . 



27r {D-l)/2 r D-l sin D-2 d 



5(r-r')5(6). 



(6.7) 



Note that the D-dimensional delta function on the right side of ( |6.6|) has been replaced 
by a cylindrically-symmetric delta function having the property that its volume in- 
tegral in D dimensional space is unity. The D-dimensional volume integral of a 
cylindrically-symmetric function f(r,8) is 



2n (D-l)/2 roc 



poo pn 

T y drr ' 1 J d9 sin D - 2 6 f{r,6). 



r(^i) 

We solve ( |6.7| ) using the method of separation of variables. Let 

g uj (r,r',6)=A(r)B(z), 
where z = cos 8. The equation satisfied by B(z) is then 



d 



d 



l-z A )—- z(D - 1)— + n(n + D - 2 



dz 2 



dz 



Biz) 



0. 



(6. 



(6.9) 



(6.10) 



where we have anticipated a convenient form for the separation constant. The equa- 
tion satisfied by A(r) is 



d 2 D - 1 d n(n + D - 2) 



dr 2 



r dr 



10 



The solution to Eq. (|6.10|) that is regular at 
polynomial f3F 



A(r) = (r ^ r'). (6.11) 
1 is the ultraspherical (Gegenbauer) 



B{z)=C^ 1+D / 2 \z) (71 = 0,1,2,3,...)- 



(6.12) 



The solution in region I to Eq. ( |6.11D that is regular at r = involves the Bessel 
function |38 



A(r)=r'- D / 2 J n _ 1+ n(\uj\r). 



(6.13) 
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In Eq. ( |6.13| ) we assume that D > 2 in order to eliminate the linearly independent 
solution A(r) = r 1 ^ l 2 Y n _ 1+ D_{\uj\r) , which is singular at r = for all n. The solution 
in region II to Eq. flb.ll ) that corresponds to an outgoing wave at r = oo involves a 
Hankel function of the first kind [BSfl 



A(r) = r 1_B/2 ^ 1+£ (|w|r). (6.14) 
Using a few properties of the ultraspherical polynomials, namely, orthonormality 
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2 1 - 2a 7rT(n + 2a) 

. i n! (n + a)r 2 (a) 

and the value of the ultraspherical polynomials at z = 1, 

r(n + 2a) 



(a ^ 0), 



d a) (r 



n!T(2a) 



(a ^ 0), 



(6.15) 



(6.16) 



as well as the duplication formula ( |2.31| ), we solve for Green's function in the two 
regions. Adding the interior and the exterior contributions, and performing the usual 
imaginary frequency rotation, we obtain the final expression for the stress P5fl: 



./ 



oo 

£ 

ra=0 
X 



[11 



1 + § )r(n + D - 2) 



2^- 1 7r^a D+1 



x-^- In ( / 

ax 



n-l+ 



„2-D 



(6.17) 



It is easy to check that this reduces to the known case at D = 1, for there the series 
truncates — only n = and 1 contribute, and we easily find 



17 2 



7T 



96a 2 



(6.18) 



which agrees with Eq. ( |2.32| ) for d = and a — * 2a. 

In general, however, although each x integral can be made finite, the sum over n 
still diverges. We can make the integrals finite by replacing the x 2 ~ D factor in the 
logarithm by simply x, which we are permitted to do if D < 1 because 



E x> Tin + a] 
V J ' = (a<0,a^-iV). 



n=0 



nl 



(6.19) 



Then the total stress on the sphere is 



oo 

rn 2 J—~< 



T(n + D 



na^ — ' n\T(D 

n=0 v 



V 



D 1 

n H 1 

2 



(6.20) 
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where the integral is 



Qu 



dx In (2xI v (x)K„(x)) . 



(6.21) 



We proceed as follows: 



• Analytically continue to D < 0, where the sum converges, although the integrals 
become complex. 

• Add and subtract the leading asymptotic behavior of the integrals. 

• Continue back to D > 0, where everything is now finite. 

We used two methods to carry out the numerical evaluations, which gave the same 
results. In the first method we use the uniform asymptotic expansions to yield 



35tt 



5657T 



VTX 7T 

^ U ~ ~2~ + 128i/ " 32768z/ 3 ' 1048576z/ 5 



+ .... 



Then using the identities 



E 



T(n + D- 2) 



71=0 



III 



E 

n=0 



T(n + D-2) 2 



for D < 2, 
for D < 0, 



(6.22) 

(6.23a) 
(6.23b) 



we obtain the following expression convergent for D < 4 (further subtractions can be 
made for higher dimensions): 



aHT(D - 1) 

oo 

+ E 

n=[N-D/2+2] 



[N-D/2+X] 

E 

n=0 

T{n + D-2) 



T(n + D- 2) 
n\ 

35tt 



vQv 



5657T 



32768z/ 2 1048576// 4 



(6.24) 



where Q u = Q u — un/2 — 7r/128z/, and the square brackets denote the largest integer 
less or equal to its argument. The infinite sums are easily evaluated in terms of 
gamma functions. 

In the second method we carry out an asymptotic expansion of the summand in 
( |6.20| ) in n, n — » oo 



T(n + D - 2) 



n! 



+ 



24D* - 176D 3 + 50AD 2 - 688D + 387 



192 



-n 



D-3 



+ . . . 



.(6.25) 
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Figure 4: Scalar Casimir stress F for < D < 5 on a spherical shell. 



The sums on n in the terms in the asymptotic expansion are carried out according to 



£ 

n=l 



D-k 



((k-D), D-k<-\. 



(6.26) 



In this way we obtain the following formula suitable for numerical evaluation: 



F — — 



i r i 



a 2 2tt 



Qd/ 



2-1 



1 ^ \r(n + D-2)v„ 1^ 
H ; r > — ; —Qu > n c k (D) 

rm-i ^ n tt^ 2^ v ; 

v 1 n=l L ' 1 



(6.27) 



where the Ck{D) are polynomials in D. 

Both methods give the same results j|9|, which are shown in Fig. ^. Note the 
following salient features: 

• Poles occur at D = 2n, n — 1, 2, 3, ... . 

• As we will see in the next plot for negative D, branch points occur at D — —2n, 
n — 0, 1, 2, 3, . . . , and the stress is complex for D < 0. 

• The stress vanishes at negative even integers, F(—2n) = 0, n = 1, 2, 3, . . . , but 
is nonzero at D = 0: F(0) = -l/2a 2 . 
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D 

Figure 5: A plot of the TM and TE Casimir stress for —2 < D < 4 on a spherical 
shell. For D < 2 (D < 0) the stress F™ (F TE ) is complex and we have plotted the 
real part. 



The case of greatest physical interest, D — 3, has a finite stress, but one 
which is much smaller than the corresponding electrodynamic one: F(3) = 
+0.0028168/a 2 . (This result was confirmed in Ref. 0.) 



The same kind of calculation can be carried out for the TM modes 50 . The TM 



modes are modes which satisfy mixed boundary conditions on the surface pq , pi 

d 



,D-2 



dr 



G(x, t; x', t') 



0. 



(6.28) 



|x|=r=a 



The results are qualitatively similar, and are shown in Fig. [|. In particular, removing 
the n = contribution from the sum of the TE and TM contributions, we recover 
the Boyer result (|4.28|) . 



7 Dielectric Ball: A Model for Sonoluminescence? 

In all the above illustrations, finiteness of the Casimir energy has been the result of 
delicate cancellation between interior and exterior modes. This cancellation occurs 
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because the speed of light is identical in the two regions. This indicates that if one 
considered, for example, a dielectric ball, the procedure given here would not yield a 
finite result, and indeed, such is the case. The same thing occurs for cavity modes, 
which are relevant for gluon (or quark) fluctuations in the bag model of hadrons — for 
a discussion of these issues see Ref . . 

A new application of such cavity modes has recently appeared, in connection 
with the suggestion by Schwinger |53| that the dynamical Casimir effect could be 



relevant for the roughly 1 million photons emitted per cycle in sonoluminescence |54 



There has been considerable controversy concerning this suggestion [|5^, |56|, |57j . The 
difficulty is that it has been impossible to extract finite Casimir energies, so that it has 
been suspected that the most divergent term, the so-called bulk energy contribution, 
could give rise to the necessary large, virtual photonic, energy needed. Dynamically, 
the difficulty has seemed to be the extremely short time scales required, which seem 
to imply superluminal velocities. 



Schwinger [|53| and Carlson et al. |56] simply use the unregulated difference be- 



tween the zero-point energy of the medium and that of the 'vacuum' within the bubble 
for their static estimates. That difference is evidently 

Ana 3 f (dk) 1 , / 1\ ,„ . 

where n is the index of refraction of the medium. This is quartically divergent, and 
with a plausible cutoff indeed can yield the required 10 6 optical photons per pulse. 
However, as we have seen, even in the simplest geometries, the naive unregulated 
zero-point energy bears no relation to the observable, finite, Casimir energy, either in 
magnitude or sign. Although Schwinger had earlier embarked on a calculation of the 
Casimir energy for a dielectric bubble, he abandoned the effort before its completion 
[|58| . However, I had considered the complementary situation, that of a dielectric ball, 
already in 1980 ||59|| . The extension to a dielectic-diamagnetic ball, with permittivity 
e' and permeability //, surrounded by a medium with permittivity e and permeability 



/i, was given in 1995 ||60|| . Using the methods given here, the following formula was 



derived for the Casimir energy: 

1 f°° °° d 



1 = 1 



where 



Si = Hx^ix) - sWMx)} 2 - e[si{x')e\{x) + sKx'hOr)] 2 , (7.3) 

where again the si , e/ are the spherical Bessel functions of imaginary argument, the 
quantity £ is 

^7 - 1 

(7.4) 



efj,> 
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the time-splitting parameter is now denoted by 8, and 



x = \y\y/en, x = \y 



IaAV- (7-5) 



It is easy to check that this result reduces to that for a perfectly conducting spherical 
shell if we set the speed of light inside and out the same, y/ejl = \/e r jj7, as well as set 
£ = 1. However, if the speed of light is different in the two regions, the result is no 
longer finite, but quartically divergent, and indeed the Schwinger result ( |7. 1| ) follows 
for that leading divergent term. 

But a divergence is not an answer. This bulk term in fact was present in the simpler 
geometry of the parallel dielectric slabs, treated in Sec. 3, where it was removed 
as a change in the volume energy of the system. In other words, the term to be 
removed is a contribution to the self energy of the material, since it renormalizes the 
phenomenologically described bulk energy of the substance. This issue is discussed 



in detail in Ref. ||61|| . Mathematically, this term is removed by subtracting off the 
uniform medium part of the Green's function in each region, just as we did in Sec. 4.1. 
We are then left with a cubically divergent energy. 

To simplify the subsequent mathematics, let us henceforward assume that the 
media are dilute and nonmagnetic, e — 1 <C 1, e' — 1 <C 1, and fj, — /i' — 1. This 
should be sufficiently good for the qualitative analysis of the situation. Then the 
subtracted energy has the divergence structure 

(7 - 6) 

Because the dimensionless time-splitting cutoff here has the structure 8 = r/a, this 
divergence has the form of a contribution to the surface tension, which therefore 
should also be removed by renormalization of a phenomenological parameter. Then, 
we are left with a finite, observable remainder. 

Previously |59], I had kept merely the leading term in the uniform asymp- 



totic expansion in estimating this finite result, believing, as with the perfectly con- 
ducting sphere, that that approximation should be quite accurate. But Brevik and 
Marachevsky calculated the next terms in that expansion and proved me wrong |62j . 
So jointly we then computed the Casimir energy exactly for a dilute dielectric sphere 
[p3| . That calculation can be done either by the point-split cutoff method sketched 
here, or directly by the zeta function trick, which swallows all the divergences. In 
the latter approach, we add and subtract the first two leading uniform asymptotic 
approximants 

i=i ^ } 

where the remainder is the difference between the exact expression and the asymptotic 



49 



one. The I sums here are evaluated according to the zeta function recipe, 



£ 

1=0 



(2 s - l)C(a) 



(7- 



Keeping only the first term in the asymptotic series gives the result that I previously 
quoted, E\ ~ — (e — l) 2 /(256a), while including the first two terms reverses the sign 
but hardly changes the magnitude |62j: E% ~ (e — l) 2 33/(8192a). Stopping here 



would give a 15% overestimate; including the remainder gives the result^ p3| 



E 



,0.004767 



(7.9) 



which is approximated to better than 2% by keeping the third term in the asymptotic 
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expansion 

What is most remarkable about this result is that it coincides with the van der 
Waals energy calculated two years earlier for this nontrivial geometry [61]. That 
is, starting from the Casimir-Polder potential [§] ( 3.37 ) we summed the pairwise 
potentials between molecules making up the media. A sensible way to regulate this 
calculation is dimensional continuation, similar to that described in Sec. 6. That is, 
we evaluate the integral 



E 



23 



vdW 



a 2 N 2 / d°r d D r'(r 2 + r' 2 - 2rr' cos 



-7/2 



(7.10) 



where 9 is the angle between r and r', by first regarding D > 7 so the integral exists. 
The integral may be done exactly in terms of gamma functions, which when continued 
7 yields ^ 



to D 



7 



E v , 



23 



dW 



15367ra 



If 



(7.11) 



This attractive result is numerically identical with the Casimir result ( |7.9| ). (Inci- 
dentally, the dilute Casimir effect for a cylinder, or equivalently, the van der Waals 
energy, vanishes [|||, |65|.) 

Thus, there is no longer any doubt that the Casimir effect is coincident with 
van der Waals forces. The fact that the latter is physically obviously proportional to 
(e — l) 2 (which had been already emphasized in Ref. ||59| ) seems convincing proof that 
the expression ( [7.1|) is incorrect. In turn, this seems to demonstrate the irrelevance 
of the Casimir effect to the light production mechanism in sonoluminescence. This 
is because nearly certainly the adiabatic approximation is valid (because the flash 
time is about 10~ n s, far longer that the characteristic optical time scale of 10~ 15 

which on the face seems to 
see Ref. |56 . 



s. (Thus the instantaneous approximation of Ref. [57 
require superluminal velocities, is not credible.) For details 



immediately after I had completed this calculation, I received a manuscript from Barton with 
the same result |64j calculated by a more elementary method using ordinary perturbation theory. 
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8 Conclusions 



In this survey of the Casimir effect, I have been very selective. Emphasis has been 
upon my own work, in part because of its greater familiarity. However, that emphasis 
also reflects my bias toward the superiority of Green's function techniques over, say, 
zeta function methods, which may reach the correct answer by defining away phys- 
ically interesting divergent contributions. The first six sections of this report treat 
unambiguously finite Casimir forces, about which there can be no argument. That is 
not to say there is not yet work to be done there, for it remains true that even the 
sign of the Casimir force is unpredictable without detailed mathematical calculation. 

Section 7 is far more speculative, in that when interior and exterior modes are 
divorced, even so mildly as through a different speed of light, the Casimir effect is 
no longer finite. Removal of that divergence requires some sort of renormalization, 
which remains controversial. Thus, it is still possible to disagree about the possible 
relevance of the Casimir effect to sonoluminescence, although I, of course, have a 
definite (negative) opinion on the matter. Recent developments in this domain of the 
Casimir effect suggest that earlier calculations of the Casimir effect in the bag model 
of hadrons be re-examined JJ7|, |45|, [52| . 



Due to lack of time and space, many other interesting topics had to be omitted. 
For example, the Casimir effect likely plays a role in the compactification of higher 
dimensional unified models It may also be relevant in 2-dimensional condensed 
matter systems p9"|. Most important are the small, but conceptually vital, inclusion 



of radiative corrections, for everything discussed here has been at the one-loop level 



70 1 . For these, and many other applications and technical developments, the reader 



is referred to the original literature, or to my forthcoming book on the subject JTT 
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